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ABSTRACT 

We present a detailed and homogeneous analysis of the ROSAT PSPC surface bright- 
ness profiles of 36 clusters of galaxies with high X-ray luminosity (Lx ~ 10 45 erg s _1 ) 
and redshifts between 0.05 and 0.44. Using recent ASCA estimates of the tempera- 
ture of the gas for most of the clusters in the sample, we apply both the deprojection 
technique and model fitting to the surface brightness profiles to constrain the gas and 
dark matter distributions under the assumption that the gas is both isothermal and 
hydrostatic. 

Applying robust estimators, we find that the gas fraction within rsoo of the clusters 
in our sample has a distribution centred on / gas (^50o) — 0A68h^ o . The gas fraction 
ranges from 0.101 to 0.245 at the 95 per cent confidence level. The values of / gas show 
highly significant variations between individual clusters, which may be explained if 
the dark matter has a significant baryonic component. Within a cluster, the average 
radial dependence of the gas mass fraction increases outward as r s , with s ~ 0.20. 
Combining these results with those of primordial nucleosynthesis calculations and the 
current estimate of Hq, the above central location implies fio,m ^ 0.56 at the 95 
per cent confidence level. This upper limit decreases to 0.34 if we take the highest 
significant estimates for / gas . 

A significant decrease in cluster gas fraction with redshift from the local value, 
/gas,o> of 0.21, found assuming ilo.m = 1, is also reduced if fio.m is low. 

Key words: galaxies: cluster: general - galaxies: fundamental parameters - inter- 
galactic medium - X-ray: galaxies - cosmology: observations - dark matter. 



1 INTRODUCTION 

The physics of the formation of clusters of galaxies depends 
upon the cosmological parameters, Qo,m and flb, that de- 
scribe the observed total matter density and its baryonic 
contribution, respectively. Assuming that clusters maintain 
the same ratio f2b /fio.m as the rest of the Universe, a mea- 
sure of the cluster baryon fraction can be compared with 
calculations from cosmic nucleosynthesis considerations of 
the abundance of the light elements (e.g. D, 3 He, 4 He, 7 Li) 
to give a direct constraint on fio, ra- 
in recent years, White et al. (1993), White & Fabian 
(1995), David, Jones & Forman (1995) and others have dis- 
cussed this issue first for the Coma cluster and then for 
samples of clusters, highlighting the necessity of a low den- 
sity Universe in order to reconcile the baryon fraction of 
the total cluster mass with the primordial fib ~ 0.05. This 
"Baryon Catastrophe" for a flat Universe persists on super- 
cluster scales; Fabian (1991) and Ettori, Fabian & White 
(1997) estimate a 15 per cent gas contribution to the total 



mass of the Shapley Supercluster, over a region of 30 Mpc 
in radius. 



The main baryonic component of the richest clusters is 
the intracluster gas, whereas galaxies contribute less than 
about 4 per cent with respect to the total gravitating mass 
(White et al. 1993, Fukugita, Hogan & Peebles 1998). Apart 
from a brief discussion in Section 5, we will not consider 
further any other contribution to the total baryon budget in 
clusters, such as baryonic dark matter or cool gas, for which 
the uncertainties are still large (cf. Fukugita et al. 1998). 
Thus, the gas fraction provides a lower limit on the total 
baryon fraction in clusters. 



The information on the gas and total mass distribu- 
tions is inferred from spectral and spatial analyses in the 
X-ray waveband. To properly know the gas density, we 
need to deproject the observed surface brightness profile, 
which is simply the projection on the sky of the (mostly) 
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bremsstrahlungjj emissivity, e oc (cluster plasma tempera- 
ture, T gas ) 1/2 x (gas density, p gas ) 2 (for T gas £ 3 x 10 7 K; 
see e.g. Rybicki & Lightman 1979), i.e.: 

f°° edr 2 

Due to the small dependence of S(b) on the temper- 
ature (in particular in the ROSAT waveband), p gas is well 
constrained from eqn. Assuming that the hydrostatic equi- 
librium holds in the cluster regions examined, we can write: 

1 rfP gas _ d(f> _ GMtotjr) (2) 
p gas dr dr r 2 

where G is the gravitational constant, and the gas pres- 
sure, P gas , is calculated through the perfect gas law, P gas = 
p ga sfcTg a s/(/im p ) (the mean molecular weight, (x, is 0.6 in 
atomic mass unit). At the present, there are two unknown 
quantities: the temperature profile (for sake of simplicity 
we assume the gas to be isothermal, but see Sect. 4.3) and 
the dark matter distribution. Fixing one of these allows us 
to solve the differential equation for the other one. In par- 
ticular, according to the different cases that we discuss in 
Sect. 3, we adopt the dark matter density profile found in 
N-body simulations and the best available estimate of the 
intracluster temperature, T gas . 

In this paper, we present the analysis of ROSAT Posi- 
tion Sensitive Proportional Counter (PSPC) surface bright- 
ness profiles of 36 clusters of galaxies, with X-ray luminosity 
greater than 10 45 erg s _1 and redshift in the range 0.05-0.44. 
These physical characteristics allow them to be well covered 
by the PSPC field of view, with the surface brightness profile 
extending to about the virial radius. 

Due to the energy-limited range of the PSPC (0.1-2.4 
keV) , we use observations of Tg as from recent published work 
on ASCA data (0.5-10 keV). 

The paper is organized as follows. In Section 2, we 
present the cluster sample and data reduction methods. In 
Section 3, we obtain constraints on the cluster dark matter 
after comparing the deprojection analysis with a straightfor- 
ward fitting approach, using both a /3-model and a gas pro- 
file obtained through the Navarro-Frenk- White dark matter 
profile. In Section 4, the value and the distribution of the gas 
fraction, / gas , are discussed. The constraints that we place 
on f2o,m using the best estimate of / gas , and the primor- 
dial nucleosynthesis results, are presented in Section 5. We 
summarize our main results in Section 6. 



2 THE SAMPLE 

We have compared samples of clusters of known X-ray lu- 
minosity (David et al. 1993, White et al. 1997, Markevitch 
1998, Allen & Fabian 1998) with the R OSA T archive in order 
to select observations of bright luminous clusters (Lx > 10 45 
erg s _1 ) at moderate redshifts (z > 0.05). This optimises 
analysis of the distribution of the gas in the outer regions of 
each cluster. 



* We adopt a MEKAL model for our analysis as described in 
Sect. 3 



Furthermore, we have not considered clusters that, al- 
though matching the selection criteria presented above, ei- 
ther have evidence of a major merger that affects both the 
determination of the centre and the hydrostatic equilib- 
rium, like Cygnus-A (Owen et al. 1997), A754 (Henriksen 
& Markevitch 1996), A2255 (Davis & White 1998), A3667 
(Rottgering et al. 1997) or are part of larger and more com- 
plex system (e.g. A85, Durret et al. 1998). 

Here we note that this sample is not complete in any 
sense: for example, if we consider the X-ray-brightest Abell- 
type clusters sample (XBACs, Ebeling et al. 1996), that is 
complete at the 80 per cent up to redshift of 0.2 for flux in 
the ROSAT hand (0.1 - 2.4 keV) larger than 5 x 10~ 12 erg 
cm~ 2 s _1 , and make the same selection done here, we find 
110 items. Of these, 28 are in common with our sample, 53 
are not available in the ROSAT archive, 19 have only HRI 
images, 8 have PSPC data but no information on the gas 
temperature and 2 are now confirmed with Lx.bol < 10 45 
erg s^ 1 . 

In Table 1, we present the list of the selected clusters 
with their basic physical parameters. The intracluster tem- 
peratures, Tgas, come from published spectral analyses. The 
clusters in our sample are very luminous hot objects for 
which only observatories with a wide X-ray energy band, 
like GINGA and ASCA, can properly measure temperatures 
using the hard tail of the X-ray spectrum. We therefore use 
ASCA measurements for all the clusters apart from A483 
and A2507 {GINGA data), A3112 (EXOSAT), A545, A2244 
(Einstein MPC), and A3888 (from its optical velocity dis- 
persion) . 

In particular, we are interested in the temperature of 
the bulk of the cluster gas, possibly not affected from the 
presence of any cooling central gas, which can lead to a 
lower (emission-weighted) temperature. Thus, in the follow- 
ing analysis, we consider gas temperatures that have mea- 
sured either excluding the core region (Markevitch 1998) or 
including a cooling flow component in the spectral fit (Allen 
& Fabian 1998). 

The error bars on the temperatures quoted in Table 1 
are symmetric la uncertainties. When the published source 
reports asymmetric errors at the 90 per cent confidence level, 
we consider the largest value and divide by 1.64 (this as- 
sumes that the errors are Gaussian). 

The PSPC images have been constructed from counts 
in the 0.5-2 keV band, where the Galactic and particle back- 
ground are minimum, after correction for instrumental and 
telemetry dead time, exclusion of times of high background 
counts. We also require the Master Veto count rate to be less 
than 170 counts s _1 (cf. guidelines for reduction of PSPC 
data in Snowden et al. 1994). Using the ROSAT Interac- 
tive Data Language (IDL) user-supplied libraries, we have 
divided the images by the respective exposure maps in the 
same energy band and corrected them for vignetting and 
exposure-time. The region of the detector support rib has 
been masked out as well all detected point sources. The clus- 
ters examined do not show evidence for a major merger, but 
may be affected by low level substructure. These clumps 
have also been masked despite their small contribution to 
the total flux. 

The surface brightness profiles have been extracted by 
estimating the X-ray centre from the exposure-corrected im- 
age after smoothing by a median filter of width of 5 pixels. 
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Mpc (H =50, q =0.5) 



Figure 1. The surface brightness profiles of the clusters in our sample in the 0.5-2 keV band, converted to the flux at the rest-frame 
energy band and corrected for the cosmological dimming by a factor (1 + z) 4 , are here compared with the Coma profile (dots connected 
by dotted line). 



The background has been calculated as an average of 
the counts s _1 arcmin -2 present between 40' and 45' in 
a sector which has no significant contamination from non- 
cluster emission (cf. column "bkg" in Table 1). 

The outer radius, i? ou t (Table 1), of the extracted and 
background-subtracted profile is defined as the maximum ra- 
dius where the signal exceeds twice the error present in that 
radial bin. This error is calculated by adding in quadrature 
the Poisson errors on the photon counts both in that bin 
and the background. 

The bin size, of at least 30 arcsec to avoid any effect of 
the Point-Spread-Function, is chosen to improve the signal- 
to-noise ratio in the outskirts of the cluster. In clusters where 
an improvement is significant, we quote in Table 1 values for 
bin sizes larger than 30 arcsec. 

Fig. |l| plots all the surface brightness profiles (in counts 
s" 1 arcmin -2 ) against radius (in Mpc). We convert from the 
angular size to the physical dimension in each cluster, using 
the following equation for angular distance 



r(Mpc) = 87.21 r(arcmm) .(3) 

We use, as cosmological parameters (Ho, go), the values (50 
km s" 1 Mpc -1 , 0.5). 

Here we note that this proper radius is inversely pro- 
portional to Ho and depends slightly upon go (a variation of 
5 per cent is observed on changing go from 0.5 to 0.01 in the 
redshift range 0.01-0.4, with the larger deviation of about 
10 per cent at the highest redshifts). We discuss further in 
Section 4 the cosmological dependence of the proper radius, 
r. 



3 DEPROJECTION AND FITTING ANALYSIS 

We analyse the surface brightness profiles to determine the 
distribution of the gas density in the clusters through the 
two usual techniques, i.e. by fitting the surface brightness 
profile and by deprojecting it. Both techniques make the 
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Table 1. ROSAT observation summary. Reference: [1] Allen & Fabian (1997); [2] Bohringer et al. (1998); [3] David et al. (1993); [4] 
Fukazawa et al. (1998); [5] Markevitch (1998); [6] from T x - cr opt relation in White, Jones & Forman (1997) and cr opt = 1307 ± 100 km 
s _1 (Girardi et al. 1997). Note: A2142 = rp800233 +wp800551n00 +wp800096 +wpl50084. iJ out is in Mpc, 'bkg' is in units of 10 -4 
count s _1 arcmin - 2 . ' This value (in units of 10 20 atoms cm -2 ) is the average of the HI column density, calculated within 1 degree and 
weighted with respect to the distance from the quoted coordinates (from the map by Dickey & Lockman 1990). 
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reasonable assumption that the observed projected cluster 
emission is due to X-ray emitting gas which is spherically 
symmetric. In particular, the former assumes a model for 
the gas density, projects it on the sky and fits it to the data 
to constrain the parameters of the model; the latter, the 
deprojection technique, makes a proper geometrical depro- 
jection of the profile and determines the gas density and 
temperature profile, assuming a dark matter distribution. 

In the next two subsections, we present these two meth- 
ods and the results obtained on the distribution of the in- 
tracluster gas. In the last subsection, the constraints on the 
dark matter in the case of the hydrostatic equilibrium are 
discussed. 



3.1 Fitting approach 

We have adopted here the following two models for the gas 
density and the projected surface brightness profile: 



(i) a single /3-model (Cavaliere & Fusco-Femiano 1976), 



Pgas = PO 



1 + 



Sb = S'o 



1 + 



2 -i 0.5-3/3 



;(4) 



(ii) a gas density profile obtained using the Navarro, 
Frenk & White (1997; hereafter NFW) dark matter pro- 
file in the hydrostatic equation where the gas is assumed 
isothermal (cf. Appendix; a first application on the Perseus 
cluster is presented in Ettori, Fabian & White 1998; a first 
theoretical discussion is found in Makino et al. 1998): 



pgas = C8o(l + x) V ^ X 



(3) 



where x = r/r B , rj — 4-KGp B r B /j,m p /(kT g!iB ) and p B — 
Pc5c(l+^) 3 f2o/^z, with (5 C equal to the characteristic density 
of the cluster and p c to the critical density (see Appendix 
in Navarro, Frenk & White 1997). 

The surface brightness profile is then obtained by numer- 
ical integration of the gas density through equation ^. 
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All the fits have been performed by a non-linear least 
squares technique (routine Curvefit in IDL). This algorithm, 
however, is sensitive to large departures for a (generally 
small) number of data points, the so-called outlier points. 
A statistical estimator that is able to properly weight these 
points is qualified as robust. Thus, we have also tested our re- 
sults with a robust estimator of the minimization of the func- 
tion (y[raw data] — y[nt])/<j[raw data], namely the downhill 
simplex method implemented in the amoeba routine (Press 
et a. 1992; IDL vers. 5.0). Generally, the agreement is good 
(with a deviation of the best estimate of the parameters of 
about 1 per cent in average) due to the large number of 
radial bins available. 

In the following analysis, we present the results obtained 
from the robust estimate of the model parameters. We show 
in Fig. |2] the count distribution of the values of the /3-model 
parameters, (r c ,f3), and for the NFW gas profile, (r s ,r/). In 
this figure, we overplot the best-fit results obtained over the 
radial range [0, Rout] and [0.2, R ou t] Mpc. The choice of the 
second radial range avoids any cooling flow (if present in the 
inner part of the cluster) on the estimates of the parameters. 
Doing this, we measure [average, dispersion] best-fit values 
of r c = [0.29, 0.19] Mpc, fi = [0.72, 0.09], and r s = [0.95, 

0. 671 Mpc, r) = [10.29, 1.55]. 

A range of tests to check whether the two sample pop- 
ulations, histograms for which are plotted in Fig. ^, have 
significantly different mean and/or variance, shows disagree- 
ment at 95 per cent confidence level of the variance ("F- 
variance test") in the populations of (3 and r\. 

With the intention to link the results of the fit using 
both the /3-model and the NFW profile, we perform a linear 
unweighted polynomial fit on the grid of parameters (Fig. |3|). 
We obtain r s = 3.17 r c and rj = 14.34 /3. These correlations 
are consistent with the best-fit results, when one model is 
fitted with the other (cf. also Makino et al. 1998). 

In the following analysis, we use the results from the 
fit with the NFW gas profile. A comparison between the 
two fits shows that the NFW gas profile provides a \ 2 lower 
than the value from the /3— model in 19 out of 36 clusters. 
Furthermore, even if the fit with the NFW gas profile is 
computationally more expensive, we show in the Appendix 
that its best-fit parameters are directly linked with the scale 
radius and normalization of the dark matter profile obtained 
from N-body simulations. 

In particular, as we discuss in Section 3.3, we are inter- 
ested in defining in each cluster the radius, rsoo, where the 
overdensity of the dark matter with respect to the average 
value is 500. Hence, after we have measured the best-fit pa- 
rameters in [0.2, Rout] Mpc, we calculate rsoo and, for only 
those clusters in which 7-500 < Rout, reiterate the fit over the 
range [0.1, 1] rsoo until the convergency on rsoo is reached, 

1. e. the estimated value of rsoo is within 0.01 Mpc. A maxi- 
mum of 3 trials is required to converge. 

The corresponding best fit values are quoted in Table 2 
and then used to describe the gas distribution through the 
deprojection of the central density. Using XSPEC (vers. 10, 
Arnaud 1996), we have converted the fitted count rates to 
the flux due to thermal X-ray emission from an optically- 
thin plasma [MEKAL code, based on the model calculations 
of Mewe and Kaastra (Kaastra 1992) with Fe L calculations 
by Liedahl (1995)], assumed isothermal at the temperature 
given in Table 1, with a fixed metallicity of 0.4 times so- 




0.5 



1.0 



1.5 



2.0 



2.5 



R out 



Figure 4. The histogram for Rout /rsoo provides an indication 
of the region 'seen' with respect to the extrapolated one. Values 
for this ratio greater than, or equal to, 1 confirm that our profiles 
can trace the mass within an overdensity in dark matter of 500. 
This happens in 23 of the 36 clusters (64 per cent). 



lar abundance, and with Galactic absorption included (cf. 
Table 1). 



3.2 Deprojection procedure 

In the deprojection technique, the count emissivity in each 
radial volume shell is calculated analytically and compared 
with the predicted counts from the emission predicted from 
an optically thin gas (described by a MEKAL model) ab- 
sorbed by intervening matter (JVh from Table 1) and con- 
volved with the response of the detector. After selecting 
boundary condition (i.e. the pressure in the outermost bin 
which allows the resulting deprojected temperature profile 
to match the observationally determined cluster tempera- 
ture), the gas temperature and density profiles are obtained, 
once a model for the dark matter distribution is defined (see 
White, Jones & Forman 1997 for a detailed discussion on 
the deprojection technique). 

The gravitational potential is described by a func- 
tional form given by the sum of two contributions: the cen- 
tral galaxy potential parametrised by a de Vaucoulers law 
(1948), assumed fixed for all the clusters and with param- 
eters -Reff = 30 kpc and velocity dispersion of 300 km s _1 
(Malumuth & Kirshner 1985), and a Navarro- Frenk- White 
potential for the general cluster. The latter potential has a 
scale parameter, core, and a velocity dispersion, ctdm, which 
is well represented, under the isothermal assumption, by the 
intracluster temperature [ctdm = (fcT gas //jm p ) - 5 ; cf. Ta- 
ble 2]. 

For some clusters with T gas larger than 10 keV, for 
which the temperature determination is not very precise 
(e.g. A2163, A2744), we fix ctdm to the optical estimate of 
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Figure 2. Histograms of the parameters from the /3-model and NFW gas profile. The histograms of the best-fit parameters in the range 
[0.2 — Rout] and [0 — -Rout] are with solid and dashed line, respectively. The median, mean and standard deviation for each distribution 
is quoted in the title of each plot. 



the velocity dispersion and require a flat temperature pro- 
file. The required values for the velocity dispersion are larger 
than the predicted one from the isothermal relation by 15 
and 50 per cent, for A2163 and A2744, respectively. 



3.3 The total gravitating mass 

From the fits on the surface brightness profile, and the de- 
projections, we obtain constraints on the gas density distri- 



bution under the assumption of an isothermal plasma in a 
spherically symmetric cluster. 

To do this, we have assumed a functional form for the 
total matter density and applied the hydrostatic equilibrium 
written in equation |^. The distribution of the total gravitat- 
ing mass is then obtained through the radial integration of 
the total matter density profile. 

We recall that our estimates of the total mass, Mtot(< 
r), are based upon (i) radial surface brightness profiles that 
reach a median i? ut of 1.93 Mpc, (ii) gas temperatures that 
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Figure 3. Correlations between the scale and slope parameters in the /3-model and NFW gas profile. 




I i i i i i i I i I 

1 

M 500 (10 15 M ) 

Figure 5. This plot shows the gas temperature vs. M500 with the 
best-fit power law functions: (dotted line) M500 = 1.67(±0.02) 
(T/10 keV) 15 10 15 M Q (x 2 = 212, with 29 degree-of-freedom) ; 
(dashed line) M 500 = 1.86(±0.04) (T/10 keV) 1 ' 93 * ' 09 10 15 M Q 
(x 2 = 184, with 28 degree-of-freedom). The clusters at high red- 
shift, with respect to the median value, are the solid symbols. 



both are corrected for any cool component in the core and 
were measured from broad-band detectors that describe well 
the shape of the thermal emission from the clusters, (iii) the 



assumption that the gas is isothermal (but see considerations 
in Sect. 4.3). 

This last assumption seems reasonable and conservative 
from recent results on simulations of gas dynamics (Evrard, 
Metzler & Navarro 1996), if one considers regions of the 
clusters where the overdensity of dark matter with respect 
to the average value is 500. Thus, we quote in the following 
analysis the value of the radius, rsoo, where this overdensity 
is reached, and M to t(< r 5 oo) = M500 using the best-fit pa- 
rameters of the NFW gas density profile and the relations 
discussed in the Appendix. 

Considering that we obtain a median value for rsoo of 
1.66 Mpc to be compared to a median R ou t of 1.93 Mpc (cf. 
Fig. ^ ), the proper gravitating mass for each cluster is here 
determined in a robust way, generally without any extrap- 
olation or application of either the r — T gaa or M to t — Tgas 
relation (Evrard et al. 1996, Hjorth, Oukbir & van Kampen 
1998), which would need a proper calibration through in- 
dependent measurements of the mass (e.g. hydrodynamics 
simulations, gravitational lensing) and replaces the peculiar- 
ity of a cluster with an average behaviour. 

For comparison with previous work, we have also inves- 
tigated the Af5oo — Tgas relation, after calculating M500 for 
the 30 clusters of our final sample (see next section). Adopt- 
ing the scaling relation M oc r T szs oc T g 3 / s 2 (l + zr 3/2 , where 
the last step implies that r oc M 1//3 (l + z)^ 1 , we obtain a 
best-fit of M 500 = 1.67(±0.02)[Tio, g as/(l + 2)l 1 - 5 x 10 15 M Q , 
where Tio, gas is in units of 10 keV (Fig. El). Leaving as 
free parameter the slope of the temperature, we measure 
M500 = 1.86(±0.04)[Tio, g as/(l + ^)] 1£,3±0 - 09 , slightly steeper 
than a power law with index of 1.5. 

At rsoo, Evrard (1997) measures from simulated clusters 
a coefficient 2. 22 ±0.32 for the above relation. The disagree- 
ment with our best-fit result is due mainly, apart from the 
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Figure 7. The gas mass fraction values for our 36 clusters are 
plotted versus the ratio / ga s/e(/gas)- Only the 30 clusters with 
/gasA(/gas) > 2 have been considered in our final sample. 

large scatter, to the steeper gas profiles predicted by the 
simulations. 



Figure 8. Bayesian probability distribution for the gas fraction 
observed in the clusters in our refined sample. This distribution 
peaks at / gas = 0.168. The arrow indicates the constraints from 
the primordial nucleosynthesis value for the low D/H case. It has a 
probability of 7.2 X 10~ 3 with respect to the plotted distribution. 



4 THE GAS FRACTION 

The standard primordial nucleosynthesis theory indicates 
the present universal ratio of baryons to photons as the free 
parameter to be constrained through the observations of the 
abundances of the light-elements. Once the microwave back- 
ground temperature is fixed, and models with three light 
neutrinos adopted, one can estimate the baryonic mass den- 
sity in units of the critical density, fib = pb/pc- Recently, 
conflicting estimates of the abundance of deuterium, D, have 
raised questions on the robustness of a general value for fi b 
(e.g. Hogan 1997). A high D abundance agrees with 4 He 
and 7 Li measurements and with the baryonic cosmic budget 
(Fukugita et al. 1998) and constrains Ob^io between 0.020 
and 0.064 (with a central value of about 0.04; e.g. Songaila, 
Wampler & Cowie 1997), on the other hand a low D/H ratio 
match Galactic chemical evolution and local measurements 
better. In this case, Q b hl is 0.076 ± 0.004 (cf. Buries & 
Tytler 1997). 

Thus, if regions that collapse to form rich clusters in 
an Einstein-de Sitter Universe retain the same value of 57b 
as the rest of the Universe, only a few per cent of cluster 
masses can be due to baryons (mostly gas in the ICM, but 
also stars in galaxies and eventual dark and cool baryons), 
in opposition with the observed 10-30 per cent (Briel et al. 
1992, White et al. 1993, White & Fabian 1995, David, Jones 
& Forman 1995, this paper). 

It is worth noting that, historically, X-ray observations 
have always shown a relatively high baryon fraction in clus- 
ters (e.g. Stewart et al. 1984); White & Frenk (1991) high- 



0.30 




0.05 0.10 0.15 0.20 
f gas (0.3 r 500 ) 



Figure 9. The gas fraction compared at different radii, i.e. at 
0.3 and 1 rsoo- This corresponds to about 0.5 and 1.7 Mpc, if we 
consider the median value of our sample of 36 clusters. The solid 
line indicates / ga s(0.3 rsoo) = /gas(1.0 rsoo). Most of the clusters 
show an excess in the gas contribution at larger radii. The clusters 
at high rcdshift, with respect to the median value, are the solid 
symbols. 
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Figure 6. Plot of the difference in standard deviation between / gas from published sources and our estimate for each cluster in our 
sample. We compare the data at the same radius and using the same technique [i.e. the deprojection for the results from Allen & Fabian 
(1998) and White et al. (1997)]. For the White et al. (1997) results, we propagate the uncertainties on the "rcprojectcd" gas temperature 
to the total mass. We include also the Coma cluster as comparison. Note: •, White et al. (1993; at r = 3 Mpc); +, A1795: Bricl & Henry 
(1996; at r = 1 Mpc), A2029: Sarazin et al. (1998; at r = 1.88 Mpc), A2163: Elbaz, Arnaud & Bohringer (1995; at r = 1.5 Mpc), A2256: 
Markevitch & Vikhlinin (1997; at r = 3 Mpc), A2390: Bohringer et al. (1998; at r = 1 Mpc); *, Allen & Fabian (1998; at r = 0.5 Mpc); 
O, White, Jones & Forman (1997, up-to-date version of White & Fabian 1995; at r = 1 Mpc); A, David, Jones & Forman (1995; at 
r = r 50o); n i Buote &: Canizares (1996, but after corrections on the total masses [Buote, priv. comm.]; at r = 1 Mpc). 



lighted the discrepancy for the Coma cluster when new 
tighter and lower constraints from nucleosynthesis were pub- 
lished by Walker et al. (1991). 

In our case, the gas fraction, / gas , is estimated by the 
ratio of the gas mass determined in the deprojection (and 
fitting) analysis to the gravitational mass estimated through 



the hydrostatic equation. The measured values of / gas at 0.5, 
1, 1.5 Mpc (where available, i.e. for radii less than Rout) and 
rsoo (from the best-fit parameters of the NFW gas density 
profile) are quoted in Table 2. 

All the values quoted at r 500 are obtained using the 
best fit parameters as described in Sect. 3.1. The 68 per 



© 0000 RAS, MNRAS 000, 000-000 



10 



Table 2. This table contains the dark matter parameters used in the deprojection analysis, the best-fit parameters from the fitting 
procedure, the value of rsoo as measured applying the NFW gas profile, and the gas fraction from the deprojection analysis calculated, 
where available (i.e. for radii less than Rout; cf. Table 1), at the bin closest to the quoted radius. In italics we quote / gas as measured at 
that radius by the fitting method. The value of / ga s(?"50o) is obtained by using the best-fit parameters with the NFW profile. Note: core, 
r c,?*s,f*500 m Mpc; (Tdm m k m s~ 1 . In the first column, we quote in parenthesis odm as calculated from T gas assuming the isothermal 
condition, t from optical analysis in Struble & Rood (1991; A2163), Smail et al. (1997; A2744). 



cluster DEPROJ FITTING / gas 

core, <t D m r s,V r c ,/3 rsoo < 0.5 Mpc < 1.0 Mpc < 1.5 Mpc < r 500 
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cent uncertainty on / gas is obtained by propagating the er- 
rors on the gas mass and total mass. The former are ob- 
tained in the deprojection analysis, by perturbing the sur- 
face brightness profile 100 times, according to the Poisson 
error on the counts in each radial bin. The uncertainty on 
the gravitating mass comes from assuming the above de- 
pendence upon the gas temperature, M to t oc Tgls, and 
propagating the relative error on the temperature itself, i.e. 
e(M tot )/M tot = 1.5e(T gas )/T gas . 

To quantify the deviation between the deprojected and 
fitted results, we calculate at R out for each cluster the quan- 
tity (/gas.DEPROJ - /gas,FiT)/e(/ g as)- On the sample of 36 
clusters, we measure a median deviation of +0.48e(/ gas ), 
with a range of —1.93, +2.81. The largest deviation is due to 
A2744, which shows a minor merger in the X-ray image and 
strong disagreement between the X-ray temperature and op- 
tical dispersion when isothermality is assumed. 



4.1 Comparison with previous work 

Our sample of clusters with high X-ray luminosity contains 
several clusters already analyzed, with the Einstein and 
ROSAT HRIs. Only 6 of them have been studied also with 
the PSPC (i.e. A401, A1795, A2029, A2163, A2256, A2390). 

We recall that the use of PSPC profiles extracted to 
about 2 Mpc (i) does not allow a good resolution in the 
inner part of the clusters, both in the deprojection analysis 
and in the fitting procedure, where we also cut the profile 
below 0.2 Mpc (or 0.1 rsoo) to avoid the contribution from 
any cooling flow, and (ii) mainly weights the outskirts in 
the fitting analysis (just a few bins are located in the core). 
These points have to be borne in mind when we compare 
our results on the gas fraction measured in the inner part 
(e.g 500 kpc) with estimates obtained through observations 
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Figure 10. The average at different radii of / gas from all 30 
clusters with the respective error. We overplot the best-fit power 
law (index of 0.20 ± 0.02; \ 2 = 0.77 for 12 degrees-of-freedom) . 



Figure 12. The gas fraction at rsoo is shown with the corre- 
sponding gravitating mass at the same radius. The clusters at 
high redshift, with respect to the median value, are the solid sym- 
bols. 
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Figure 11. This plot shows the dependence of / gas on the 
redshift with the best-fit result of 0.215(-0.019, +0.020) (1 + 
z )-l.75(-0.65,+0.6B) ^2 = 1Q2) w ith 1 degrees-of-freedom) . The 
rr-position of the bins is a biweight location. The y-position is the 
biweighted value with the respective error. When a different cos- 
mology is considered into the dependence of / gas upon rf a ^ 5 g , we 
obtain the values represented with diamond (Qo >m = 0.2, = 0) 
and triangle (f2o, m = 0.2, SIa = 1 — ^o,m)- The dotted line rep- 
resents the central value of / gas from the robust analysis on 30 
clusters. 
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1 — fio.m- 
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which are more efficient for that purpose (cf. the results from 
HR1 observations in Allen & Fabian 1998) . 

In Fig. ^, we plot the difference in standard deviations 
between previous estimates and our results from the fitting 
procedure. We generally underestimate the gas fraction in 
the core of cooling flow clusters (e.g. A478, A1689, A1835, 
A2029, A2142) when compared to the Allen & Fabian (1998) 
results. 

In A1795, the effect of the cooling flows produces a re- 
markable mismatch between our estimates and the results of 
David et al. (1995; about 3.7a) and Allen & Fabian (1998; 
~ 5.5<t). On the other hand, the conclusion from Briel & 
Henry (1996) shows an opposite trend with respect to these 
and to our result from the fitting procedure, while White et 
al. (1997) are in good agreement with our value at 1 Mpc. 



4.2 The distribution of / gas and its dependence 
upon redshift 

We select from our sample the clusters where the value of 
/gas(< ?"5oq) is twice its uncertainty, i.e. we apply the selec- 
tion criterion that / ga s/e(/gas) > 2 (Fig. ^). This leaves 30 
clusters (A545, A586, A2244, A3112, IRAS09104, MS1358, 
excluded) as our definitive sample. As discussed in Sect. 2, 
this sample is not complete in any sense. In the following 
analysis, we adopt the null hypothesis that (i) the gas frac- 
tion is constant with redshift and total mass, (ii) our high 
X-ray luminosity clusters are a homogeneous sample where 
no contamination is expected by, say, cooling flows. (Note 
that the gas temperatures are obtained from spectral anal- 
yses which exclude the cooling flow regions.) 

The weighted mean of the values of / gas (r5oo) is 
0.176(±0.003), with a variance of 0.037. 

The mean, however, is statistically a poor estimator, 
being sensitive to (i) the replacement of even a small part of 
the data with different values, (ii) the assumed underlying 
population, (iii) the increase in the number of data to get 
better information. In this sense, the median is certainly a 
better estimator. But, an even better (i.e. more resistant, ro- 
bust and efficient) estimator is the biweight location (Beers, 
Flynn & Gebhardt 1990). Using this estimator, we calculate 
/ g as(< rsoo) = 0.171 ±0.035. Moreover, we can use the clus- 
ters where a significant cooling flow is present (cf. Table 1 
in Allen & Fabian 1998) to estimate the gas mass fraction 
characteristic of more relaxed systems. We select 13 clusters 
(A478, A1068, A1413, A1689, A1795, A1835, A2029, A2142, 
A2204, A2390, MS2137, PKS0745, Zw3146) and measure a 
biweight of 0.168 ± 0.030, that is consistent with the distri- 
bution in the whole sample. 

Both the weighted mean and our estimate of the bi- 
weight location are completely consistent with the biweight 
estimate of the gas fraction, at rsoo, of 0.170 ±0.008 quoted 
by Evrard (1997) and calculated from two different pub- 
lished samples (White & Fabian 1995, David et al. 1995) 
after revision of the total mass of the David et al. clusters. 

We now verify if the independent measurements of / gas 
are compatible in a significant way. Applying the x 2 test to 
our sample of / gas , we find that these are not compatible at 
the > 99.9 per cent confidence level when the weighted mean 
is adopted as the representative value. In fact, the deviations 
between our values of / gas (rsoo) can be up to a factor of 3 (cf. 



A483 and A2142; see also conclusions on A1060 and AWM7 
in Loewenstein & Mushotzky 1996). 

Carrying over our robust approach to these statistical 
issues, we investigate the compatibility between each value 
of / gas and its average representative estimate, following the 
considerations in Press (1996) on the Bayesian combination 
of apparently incompatible measurements of the underlying 
quantity (in our case, / gas ). This method, which weighs a 
weighted sum of 'good' and 'bad' Gaussians attributed to 
each measurement with an a priori probability that the ex- 
periment is 'correct', provides a distribution of probability 
for the value under investigation and a judgment on the 
goodness of each experiment. Applying it, we find that all 
the measurements, with a significantly low probability of 
less than 0.01 per cent that they are all wrong, are con- 
sistent with a distribution, roughly symmetric, peaked at 
0.168 and with a dispersion of of [-0.030, +0.036] (Fig. §). 
This distribution is completely consistent with the biweight 
location and scale. Hereafter, we consider this value, with 
the respective "dispersion" around the "mean", as the rep- 
resentative estimate of / gas in our sample. In Fig. [| we also 
indicate the estimated value of fib from low D abundance. 
When we compare this value with the calculated probabil- 
ity distribution of / gas (rsoo), we locate it on the wing of the 
distribution as very unlikely (probability of 0.7 per cent). 

We investigate now two other issues related to the dis- 
tribution of gas in clusters: (i) the dependence of the gas 
fraction on the radius within each cluster, and (ii) its con- 
stancy with the cosmological time. These issues are strictly 
related to the formation and evolution of clusters of galaxies 
in the present cosmological scenario. 

The observed structure in the Universe results from the 
evolution of gravitational instability. On scales between 10 13 
and 1O 15 M0, gravity is assumed to be the only force driv- 
ing the formation of groups, poor clusters and clusters of 
galaxies. In this scenario, where the evolution is an entirely 
self-similar process, the gas and the total mass should have 
the same distribution. 

This is not what is generally observed. In clusters, data 
analysis (e.g. White & Fabian 1995 and David, Jones & For- 
man 1995) and hydrodynamics simulations (Evrard 1997) 
show an increase of / gas with radius. A way to explain the 
discrepancy between the distribution of gas and the underly- 
ing dark matter is to take into account physical phenomena 
able to redistribute the energy in the cluster, like galactic 
winds, ram-pressure stripping and heat input by supernovae 
type II (Evrard 1997, Metzer & Evrard 1998, Cavaliere, 
Menci & Tozzi 1998, Wu, Fabian & Nulsen 1998). In partic- 
ular, the role of galactic winds is now well-studied in raising 
the gas entropy from the value achieved after gravitational 
collapse and flattening the gas distribution inside the clus- 
ter. This should affect both the global value of / gas and its 
radial dependence. 

In our case, which selects high-luminosity clusters, the 
median value of T gas of 9 keV should save our local / gaa esti- 
mates from dropping significantly with respect to the global 
value. Feedback from galaxy formation affects the low- 
temperature clusters most strongly (e.g. Metzer & Evrard 
1998), where the total intracluster thermal energy becomes 
comparable to the energy input from feedback. 

On the other hand, we can investigate the redistribu- 
tion of cluster energy, describing the radial variation of the 
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Table 3. The biweight values of / gas (with the respective boot- 
strap error) for different subsamples selected, first, in redshift 
with respect to the median value of 0.1680, and then in M500 
(the median values for the low and high redshift sample are 0.97 
and 1.08 x10 15 Mq, respectively). In parenthesis, the number of 
values for each subsample is quoted. 



low z high 2 

low M500 0.189 ± 0.020 (7) 0.167 ± 0.012 (7) 
high M 500 0.192 ± 0.017 (8) 0.151 ± 0.012 (8) 



gas fraction. In Fig. O, we show the general increase in / gas 
when calculated at 0.3 and 1.0 rsoo. To quantify this trend, 
we collect in 14 radial bins equally spaced between [0.3, 1.0] 
rsoo the averaged gas fraction (with an error on the mean ob- 
tained from the propagation of the individual errors) , mak- 
ing a composite f gas profile. Then, fitting a power law to 
the 14 bins (Fig. [ji]), we obtain / gas (r) oc (r/r 5 oo) - 20±0 '° 2 . 
This result is slightly steeper than the estimate of about 
0.13-0.17 from simulations (cf. Evrard 1997). 

Moreover, in the two clusters (A644, A1651) for which 
the radius where a dark matter overdensity of 200 is sur- 
veyed (i.e., r2oo < -Rout), we observe that the gas fraction 
increases by about a further 15 per cent from rsoo to r2oo- 

The evidence of a positive gradient underlines a lower 
concentration of the gas with respect to the dark mass, that 
is consequence of both our assumption on the dark matter 
profile (at r ~ rsoo, Pdm oc r~ 2 ' 4 ; cf. N-body simulations 
results, for example Thomas et al. 1998) and our best-fit re- 
sults on the slope of the gas density (i.e. p gas oc r -3 ' 3 oc r -2 ' 2 , 
that implies a dependence of A/ gas on r 0,8 ). In other words, 
converting the dependence of the gas fraction upon radius to 
a dependence on the different mass components, we conclude 
that/ gas (r) oc [Af gas (r)/M gas (rsoo)] 1/4 oc [Af tot (r)/A/soo] 1/3 . 

Another argument against simple self-similar cluster 
evolution would be a variation of / gas with redshift. 

First, the use of Spearman's rank-order correlation 
shows that / gas is correlated with redshift at a confidence 
level > 99 per cent. We investigate this further by divid- 
ing the sample according to the median value in redshift 
(z—0. 1680) and calculating the biweight location (and re- 
spective bootstrap error) of the / gas (rsoo) values. We ob- 
tain 0.192(±0.012) and 0.159(±0.009), at low and high red- 
shifts, respectively. The evidence of a higher local gas frac- 
tion is also clear from Fig. [n], where we sample the data 
in 3 redshift bins of 10 elements each. Each bin is repre- 
sented with the biweight location, and the respective error. 
We fit then a power law, / gas ,o X (1 + z)~ a , and measure 
/gas.o = 0.215(-0.019, +0.020) and a = 1.75(-0.65, +0.65), 
where the errors are la deviation. Again, the slight negative 
evolution in the gas fraction appears significant at ~ 2.7<r 
level (estimated from the value of a and its error). 

We have also looked for trends of / gaa with the total 
mass, A/500, (Fig. |l2| ) considering the biweight gas fraction 
values calculated in subsamples that were selected according 
to the median of the redshifts and of the masses (Table 3). In 
particular, we search for the minimum in the \ 2 distribution, 
when the data in Table 3 are compared with the model 

/ g a S (r 5 oo) = const x (1 + z)~ ai . (6) 

We obtain a minimum \ 2 of 0.39 (one d.o.f.) with the best-fit 



parameters: Qo = 0.15, cti = 1.24. In Fig. |l^, we plot the 68 
and 95 per cent errors contour for the two slope parameters. 
The dependence on the redshift is well in agreement with the 
previous results plotted in Fig. [ll]. On the other hand, there 
is no evidence for any statistically significant dependence of 
/gas (r 5 oo) on A/500 • 

Allen & Fabian (1998) show a decrease in the gas frac- 
tion of clusters at higher temperature (cf. their Fig. 3). If 
we replace A/500 in equation |^ with the intracluster temper- 
ature, we find that the slope of T gas is within the range [-1.4 
,0.4] at the 95 per cent confidence level. But putting Qi = 0, 
ao has to be larger than (95 per cent c.L), in agreement 
with Allen & Fabian results. 

We also note that, from the scaling low between mass 
and temperature (M oc Tg[ 2 ), any apparent dependence of 
/gas upon T gaa becomes weaker by a factor 2/3 when applied 
to mass. 

This scenario implies that the gas component in the 
X-ray highly-luminous systems considered here is almost in- 
dependent of the mass and only slightly on the temperature, 
once the redshift dependence is taken into account. 

Any apparent decrease with redshift of the gas fraction, 
however, has been recently questioned. Following an origi- 
nal idea of Sasaki (1996), Cooray (1998) and Danos & Pen 
(1998; see also Rines et al. 1998) have shown the angular 
distance— redshift relation, d ang (z, qo), that we write in equa- 
tion |^, to be the major factor responsible for the apparent 

negative evolution. From the definition of the gas fraction, 

3/2 

it holds that / gas oc d ang and, consequently, / gas tends to be 
lower at higher redshifts in a high density universe. 

We consider the changes in / gaa for cosmologically dif- 
ferent scenarios in Fig. Fill. We lower f2o, m to 0.2, apply- 
ing both an open universe and a flat one with a cosmo- 
logical constant Qa = 1 — ^o,m [cf. eqn. (25) in Carroll, 
Press & Turner 1992]. In the latter case the change ap- 
pears more significant in flattening / gas to a non-evolving 
value. Fitting a power-law to these values, we obtain / gas = 
0.210(-0.020, +0.022) x (1 + 2 )-o-6i(-o.66,+o.64) with a re _ 

duced x 2 of 1.3. All the errors are one standard deviation. 

Also when we apply equation [] to the case [f2o,m, Oa] = 
[0.2, 0.8], that shows the better agreement with the constant 
/ g as assumption, we observe a significant flattening in the 
redshift dependence. The best-fit parameters (\ 2 = 0.22 for 
1 d.o.f.) are now (ao, Qi) ~ (— 0.1, 0.5), with a range for the 
Qo values of [—0.5,0.7] at the 95 per cent confidence level 
(dotted contours in Fig. [l3| ). 

Thus, a low density Universe not only matches the ob- 
served gas fraction with the baryonic amount provided dur- 
ing the primordial nucleosynthesis, but also flattens to a con- 
stant (in look-back time) gas mass fraction of about 0.196 
(biweight estimate). 

4.3 The dependence of / gas on the temperature 
profile 

In this Section, we remark on the role of a negative gas 
temperature gradient (instead of an isothermal profile) in 
estimating the cluster gas fraction. Markevitch et al. (1998) 
have recently claimed, from the analysis of a sample of clus- 
ters observed from the X-ray satellite ASCA, that a nega- 
tive temperature gradient is generally present and well rep- 
resented by a polytropic function with index, 7, of 1.2-1.3 
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Figure 14. The gas fraction relative to the isothermal case, af- 
ter correction for the gradient in the temperature parametrised 
by the polytropic equation (see Appendix). When we assume the 
isothermal gas temperature equal to the maximum value of the 
temperatures allowed from the polytropic profile (i.e. the central 
value, To), we observe an increase of the polytropic / gas esti- 
mate up to 60 per cent with respect to the isothermal / gas value 
at r ~ 6r c (solid line). Using an emission- weighted average of 
the temperatures spanned from the polytropic profile as reference 
value for the isothermal profile (instead of the central value), the 
correction on / gas is by about 17 per cent at the same radius 
(dashed line). 



between 1 and 6 X-ray core radii. They have assumed a (3- 
model with [/3, r c ] of [0.67, 0.3 Mpc] to constrain the gas 
density, and then used the polytropic equation to fit the 
temperature profiles. 

In the Appendix, we derive the equations necessary to 
describe the polytropic state of the intracluster plasma. Us- 
ing the results of Markevitch and collabo rato rs, we make a 

of about 1.2 - 



correction to the total mass via equation A9 



0.6 over the radial range [1, 6]r c on the total mass. It yields 
a larger gas fraction with respect to the isothermal case, 
above 2 r c . At r ~ 6 r c ~ rsoo, the polytropic value of / gas 
can be larger by a factor of 1.6 than the isothermal estimate 
(Fig. 0). 



We have also fitted eqn. 



A7-A5 



to the profiles of the 34 



clusters in our sample with enough data points, to obtain 
7 = 1.25 ± 0.13. None of these polytropic fits, however, is 
better than the isothermal one at 95 per cent confidence 
level using the F-test. 

We note that steep electron temperature profiles need 
not be representative of the state of the gas. As shown by, 
e.g., Ettori & Fabian (1997) and Takizawa (1998), the effi- 
ciency of the Coulomb collisions between ions and electrons 
has to be considered when a such negative gradient is ob- 
served in clusters, even if they underwent an ancient merger. 
Thus, when the equipartition time is comparable to the age 
of the last, large, merging event, it is also true that the 



mean gas temperature profile is generally flatter than the 
(emission-weighted) electron temperature profile estimated 
from X-ray observations. 

In conclusion, the presence of a non-isothermal temper- 
ature profile tends to increase the gas fraction value at rsoo, 
more or less significantly depending on our underestimate 
of the drop of the temperature in the outskirts of each clus- 
ter with respect to the assumed representative isothermal 
value. Whether this is a systematic or random uncertainty 
in our results depends upon whether all clusters have similar 
profiles or not up to rsoo- Current data are unable to clarify 
this important issue. 



5 fi , M FROM THE GAS FRACTION 

The range of observed cluster gas fractions can now be 
compared with the results of primordial nucleosynthesis 
calculations. Assuming that the current best estimate for 
H = 73 ± 14 km s" 1 Mpc" 1 (Freedman et al. 1998), and 
the constraint 



fit),] 



fib 

Jga 



Ho 



50 km s" 1 Mpc" 1 



(7) 



we obtain fi , m < 0.20(-0.13, +0.13), adopting fi b = 0.04 
which corresponds to the high deuterium abundance esti- 
mate, and fi , m < 0.37(-0.16, +0.09), using fi b = 0.076. 
The error bars come from the propagation of the accepted 
range for the given fib (see beginning of Sect. 4). From the 
latter estimate, we can put an upper limit on the density 
of the Universe as inferred from the measured gas fraction 
in our sample of highly-luminous clusters: fio, m has to be 
less than 0.56 at the 95 per cent confidence level. We note 
that our estimate for fio, m is conservative due to our use 
of an 'average' value for / gas . If we adopt the 95 per cent 
lower limit to the highest reliable values in our sample of 
gas fraction (considering 10 per cent of the sample, that of 
A2142, A2319, A401, A2029) to constrain the cosmological 
parameter, we require that fio, m < 0.34. 

Using the conservative constraint, it is straightforward 
to require fiA = 1 — fio.m > 0.44 (95 per cent confidence 
level), still marginally consistent with the recent results on 
the magnitude-redshift relation for the type la supernovae 
(Perlmutter et al. 1998). 

We have neglected so far any other baryonic contribu- 
tion, like galaxies and baryonic dark matter. The luminous 
mass in galaxies is about a fifth of the gas fraction, i.e. 



g — /gal//g; 



0.2hlj (White et al. 1993, Fukugita et al. 



1998). The cluster baryonic dark matter is suggested to be in 
the form of cold clouds or low mass stars and brown dwarfs 
deposited by cooling flows (Thomas & Fabian 1990, Fabian 
1994). Also on galactic scale, low mass dark objects can be 
deposited by a cooling flow in the halo during first collapses 
of a protogalaxies (Nulsen & Fabian 1997). 

The last equation can be then rewritten considering 
these other baryonic contributions. Defining 6 as the ratio of 
this baryonic dark matter to the gas fraction, and assuming 
it is independent of Ho, we can write: 



fio,m — 



"50 



fib 



/ gas l + (b + g)h}J 



(8) 
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Figure 15. Expectation values for f2o, m when a baryonic contri- 
bution is considered in addition to intracluster gas and stars in 
galaxies [cf. equation B. The case b = corresponds to the Oo,m 
value estimated considering only / gas and assuming a low value 
for D/H. 

For different b values, using the best estimates from our 
previous analysis and propagating the errors as usual, we 
plot in Fig. ha the estimated fio,m with la uncertainty. For 
example, recent estimate on the universal primordial baryon 
fraction from MACHO results (Steigman & Tkachev 1998), 
assumed baryonic and in the halo of the Galaxy, give us a 
value of b ~ 0.6, with a consequent fio,m of about 0.15. 

Here we note that the contribution from baryonic dark 
matter can be about 3 times the luminous matter in normal 
stars, as predicted from standard model of gas processes in 
galactic formation (Nulsen & Fabian 1997). This can par- 
tially explain the large scatter in the gas fraction distribu- 
tion (for example, cf. Fig. |l2| ), suggesting that most of the 
baryons in clusters with lower values of / gas are in the form 
of dark matter. 



(i) At the radius where the cluster overdensity is 500, the 
measured values of / gas have a biweight location and scale 
of 0.171 ± 0.035. When Bayesian statistics are adopted, the 
probability distribution is peaked at 0.168 and has a 95 per 
cent range of [0.101, 0.245] (Fig. |). 

With respect to this distribution, the highest estimate of 
fib (corresponding to a low deuterium abundance) has a 
probability of 7.2 x 10 -3 in an Einstein-de Sitter Universe 

(fi0,m = 1). 

When we consider the more relaxed clusters with a central 
cooling flow, we measure 0.168 ± 0.030. 

If we consider the cosmological correction in the angular 
size-redshift relation for a flat Universe with fiA = 0.8, we 
measure the biweight values of 0.196 ± 0.035. 

(ii) Many of the individual estimates of / gas are inconsis- 
tent with the average value or with each other, confirming 
that there are real differences in the measured gas fractions 
in clusters (Fig. (?], |J |l^). It remains possible that some of 
this spread is due to a breakdown in our assumptions. There 
may, for example, be a merger taking place along the line of 
sight in some clusters. If the gas is actually more extended 
along the line-of-sight then we have underestimated the gas 
fraction, conversely if it is foreshortened then the gas frac- 
tion is overestimated. 

(iii) The average dependence upon the radius within a 
cluster is r" , with s ~ 0.20 (Fig. llOj) . This follows from our 
best-fit results on the slope of the gas density, p gas oc r~ 2 ' 2 , 
and our assumption of the Navarro-Frenk- White functional 
form for the dark matter distribution. In general, we can 
write / gas oc r _2,2_SDM , where sdm is the slope of the dark 
matter profile. 

(iv) In our sample of high- luminosity clusters, there is 
a highly significant correlation between / ga s(>"50o) and red- 
shift, that we quantify fitting the function / gas ,o x (1 + z)~ a : 
/ gas ,o = 0.215lg;° 2 9> a = 1.75 ±0.65 (fi , m = 1). When 
we take it into account, there is then only a mildly signif- 
icant dependence of /gas(f"50o) upon the temperature (and, 
even weaker, upon the mass). The dependence upon red- 
shift weakens (i.e. the significance of any decrease in / gas 
is reduced) in a low matter density universe, in particular 
if a cosmological constant is present. The normalization of 
the power law, i.e. the extrapolated local value of the gas 
fraction, / gas ,o, remains stable at 0.21. 



6 SUMMARY AND CONCLUSIONS 

In our sample of 30 hot (median T gas = 9 keV) and in- 
termediate distance (median z = 0.1680) clusters of galax- 
ies, we have applied both the deprojection technique and 
the fitting approach, using both (a) the /3-model and (b) a 
gas profile obtained from the Navarro-Frenk- White poten- 
tial. For (a) we find that the median best-fit parameters are 
[/3,7-c] = [0.71,0.24 Mpc]; for (b), [rj, r 3 ] = [10.20,0.75 Mpc]. 

The two techniques produce the same values of / gas 
within a median deviation of less than la. 

Assuming an isothermal profile for the gas temperature, 
with the caveats discussed in Sect. 4.3 on the role played 
by a decrease of the temperature outward, we obtain the 
following results on the gas fraction, / ga3 : 



(v) Adopting both a low and a high abundance for deu- 
terium and requiring that any physical source in redistribut- 
ing the energy within the cluster cannot produce large vari- 
ations in the value of / gas quoted above [cf. (i)], we constrain 
fio,m < 0.56 at the 95 per cent confidence level. If we take 
the highest significant estimates of / gas in our cluster sam- 
ple, we find that fio.m < 0.34. 

(vi) Future X-ray missions can constrain better the tem- 
perature gradient in the electron (and electron+ion) popu- 
lation. If the indication of a negative gradient in T gas is con- 
firmed, / gas rises and fio, m drops by a factor that can be, 
at maximum, 1.6. Such a large value of / gas could become a 
real concern for the physics of the formation and evolution 
of clusters of galaxies in present cosmological scenarios. 
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APPENDIX A: ANALYTIC MODELS FOR THE 
STATE OF THE INTRACLUSTER GAS 

The Navarro, Frenk & White (NFW) density profile, 
Pnfw = p s a; _1 (l + x)~ 2 , generates a gravitational poten- 
tial of the form: 



Rodd H.J., 1991, ApJS, 77, 363 
1998, ApJ, 509, 579 
, Fabian A.C., 1990, MNRAS, 246, 
et al., 1998, MNRAS, 296, 1061 

, Schramm D.N., Olive K.A. 



— = 4nGr s p : 
dr 



ln(l + x) - x/(l + x) 



(Al) 



where, defining r s as a scale radius, x = r/r s and p s = 
/9 c 5c(l+z) 3 ^o/^z, with 8 C equal to the characteristic density 
of the cluster and p c to the critical density (see Appendix 
in Navarro, Frenk & White 1997). 

Putting it in the hydrostatic equation, and assuming 
the isothermality for the intracluster gas, T gas (r) = To, 



1 rfpgas 

p g as dr 
we obtain 



(imp 



(A2) 
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-- Po e-\l + xY /x =a (l + x)" /l , (A3) 
where 

_ 4-KGp s r*nm p _ 1.55 c (l + z) 3 (Cl /Cl z )Hor*iu,m p 

v - w - Wo • (A4) 

(we define, as usual, G as the Gravitational constant, fio 
and Q,z as the cosmological parameters at the present time 
and at redshift z, and Ho as the present-time Hubble con- 
stant). The gas density profile is then integrated and fitted 
to the observed surface brightness profile to constrain the 
free parameters po,r a and r\. Hence, given To and the best- 
fit parameters r B and rj, 5 C — 8/3 x c 3 /[ln(l + c) — c/(l + c)] 
and the concentration parameter c can be estimated from 



eqn. <\4. Here, 8 is the mean interior overdensity that de- 
fines a cluster over the background. It is generally assumed 
equal to about 200, that is the density contrast of a virialized 
object in the nonlinear regime of spherical collapse. 

The radius at which this mean overdensity is reached, 
rj, is then defined rj = c r s (Navarro, Frenk & White 
1995, 1997). NFW (1997) show that for a cluster iden- 
tified at redshift z, the cosmological dependence of rj is 
(fio/n z )- 1/3 (l + z)- 1 = (1 + n z)- 1/3 (l + z)~ 2/3 . This de- 
creases with the redshift and increases slightly with the low- 
ering of the density parameter: for example, at redshift 0.4, 
rj is lower than the local estimate by about 29 and 23 per 
cent for Q,q equal to 1 and 0.3, respectively. 

We can also generalize to an intracluster gas with a 
polytropic equation of state 



Po 



r gas (r) 

To 



(A5) 



the consequences of the hydrostatic equilibrium with a 
clust er p otential described by the NFW functional form 



(eqn. Al) 



d(fcT gas ) _ 7- 



pm p - 

dr 7 

that integrated over the radial range of interest provides the 
equation: 



To 



7- 1 
1 + J7- 



ln(l + x) 



(A7) 



where T gas (0) = T and 77 = r?(T ). 

The gas density profile is then obtained using equa- 
tion 



A5 



The presence of a temperature gradient also affects the 
estimate of the total gravitating mass, Mtot, for which the 
gravitational potential is unkwown: 

Mtot = 7 (e^Y~ 2 * ( Est). (A8) 

p,m p G \ po J dr y p J 

The case 7=1 provides the usual equation to estimate 
the total gravitating mass when the gas is isothermal. In 
particular, the ratio between the polytropic and isothermal 
mass estimates will be given by 

MP°* _Jp ga »Y~ 1 _ 7 Z-W. (A9 ) 



M 1S ° \ Po J To 

This equation has been used to estimate the correction 
on / gas for the presence of a temperature gradient as plotted 
in Fig. 14. 
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